Interaction-Enhanced Coherence Between Two-Dimensional Dirac Layers 
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We estimate the strength of interaction-enhanced coherence between two graphene or topological 
insulator surface-state layers by solving imaginary-axis gap equations in the random phase approx- 
imation. Using a self-consistent treatment of dynamic screening of Coulomb interactions in the 
gapped phase, we show that the excitonic gap can reach values on the order of the Fermi energy at 
strong interactions. The gap is discontinuous as a function of interlayer separation and effective fine 
structure constant, revealing a first order phase transition between effectively incoherent and inter- 
layer coherent phases. To achieve the regime of strong coherence the interlayer separation must be 
smaller than the Fermi wavelength, and the extrinsic screening of the medium embedding the Dirac 
layers must be negligible. In the case of a graphene double-layer we comment on the supportive role 
of the remote 7r-bands neglected in the two-band Dirac model. 

PACS numbers: 71.35.-y, 71.35.Lk, 73.22.Pr, 73.22. Gk, 68.65.Pq 



I. INTRODUCTION 

States in which quantum coherence is spontaneously 
established between electrons in two different layers are 
known as bilayer exciton condensates. This type of or- 
dered electronic state^ has so far been realized only in 
GaAs/AlGaAs double quantum wells, only at low tem- 
peratures, and even then only in the presence of a mag- 
netic fielcP. When the separate layers of the bilayer con- 
densate are contacted independently, these systems ex- 
hibit novel low bias-voltage collective transport anoma- 
lies^lrelated to interaction-enhanced interlayer coherence. 
These anomalies could enable low-power logic devices 
based on new physical principles, for example the BiS- 
FETBj device. It is therefore of interest to search for 
bilayer systems in which there is a potential for sponta- 
neous or strongly enhanced coherence in the absence of a 
magnetic field and at as high a temperature as possible. 

In the absence of a magnetic field, bilayer exciton con- 
densates are likely to form when the electron and hole 
Fermi surfaces in the two layers are nested. For this 
reason the prospects for realizing bilayer exciton conden- 
sates in the absence of a field have improved recently with 
the discovery of two new classes of Dirac two-dimensional 
electron systems, one based on graphene sheets^ and the 
other based on topological insulators^. In Dirac two- 
dimensional electron systems, electrons and holes have 
identical Fermi surfaces at equal carrier densities because 
of the linear band crossing at the Dirac point. The in- 
dependent control over the carrier type and density on 
each layer, easily realized with gate voltages, makes these 
materials attractive candidates to observe interaction in- 
duced interlayer coherence. 

The temperature at which the condensation sets in is 
difficult to estimate theoretically because, among other 
issues, it is sensitive to dynamic screening of the electron- 
electron interaction. This problem has been studied 
previously using a variety of approaches and approx- 
imations^mu. Xo our knowledge only one study has 



accounted for the dynamic nature of screening^, and 
none have accounted simultaneously for the reduction of 



screening^ that accompanies the appearance of an inter- 
layer coherence gap. 

In the present paper we show that both the dynamic 
nature of screening and its reduction in the coherent 
phase play an essential role in determining the properties 
of the exciton condensate in double-layer Dirac systems. 
Our findings suggest that spontaneous coherence is pos- 
sible at relatively high temperatures in suitable double- 
layer graphene systems. (We refer to a system with two 
weakly hybridized graphene layers as a graphene double- 
layer in order to distinguish it from a Bcrnal stacked 
two-layer graphene system, which is often referred to as 
bilayer graphene.) 

We use a simple approximation that employs a two- 
band Dirac model for the electron and hole systems and 
accounts for dynamical screening in a random phase ap- 
proximation. In single-layer graphene systems the ran- 
dom phase approximation provides an accurate descrip- 
tion^ of photoemission experiments^ and in particular 
correctly accounts for the plasmaron features seen near 
the Dirac point of systems with finite carrier densities. It 
can also be shown theoretically to apply in double-layer 
graphene systems due their large number of fermion fla- 
vors^. Here we use this approach to assess the require- 
ments for spontaneous or enhanced coherence in both 
graphene double-layer and topological insulator systems. 
We conclude that strongly enhanced coherence at room 
temperature is possible in the graphene case, provided 
that high carrier densities can still be realized at very 
small interlayer separations and the extrinsic dielectric 
screening from the substrate can largely be eliminated. 
Strongly enhanced coherence at room temperature in 
topological insulators, on the other hand, may require 
a search for new materials with larger bulk gaps. 

Our paper is organized as follows. In Sec. [TTJ we de- 
scribe the two-band double-layer model that we use for 
concrete calculations. In Sec. |IH|we discuss the proper- 
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FIG. 1: (color online) (a) Illustration of a double-layer sys- 
tem. The two layers represent either two sheets of graphene 
separated by a thin dielectric, or the two faces of a topolog- 
ical insulator thin film. The layers have carriers of opposite 
charge and equal densities, and are in thermodynamic equi- 
librium. Dirac cones are displaced by the presence of carriers 
in each layer (b), there is an interval of Fermi energy with 
only two Fermi surfaces, one from the conduction band of the 
high density layer and one from the valence band of the low 
density layer. For a particular value of Fermi energy (c) the 
electron and hole carrier densities are equal so that the two 
layers have a common Fermi surface. Under these circum- 
stances the system is unstable to a broken symmetry state 
with spontaneous interlayer coherence and a gap (d) in the 
quasiparticle spectrum. 



ties of screening in a coherent bilayer system. In Sec. |IV| 
we analyze the interaction-enhanced coherence in a con- 
stant gap approximation, which does not take into ac- 
count the wavevector and frequency dependence of the 
anomalous (electron-hole pairing) self-energy. In Sec. M 
we determine the spontaneous gap including its wave- 



vector dependence. In Sec. [VI] we discuss the role of the 
remote bands that do not intersect the Fermi surface in 
establishing the interlayer coherence. Finally, in Sec. |VII| 
we summarize our results and present our conclusions. 



II. DOUBLE-LAYER MODEL 

We consider a double-layer system with a layer sep- 
aration d. Each layer hosts a two-dimensional electron 



gas that is described by a Dirac model with velocity vd- 
We assume that the Fermi levels ep have been adjusted 
so that the magnitude of the carrier density is the same 
in both layers as illustrated in Fig. [T] The conduction 
band of the electron-doped layer and the valence band 
of the hole-doped layer therefore share a common Fermi 
radius kp, whereas the remote bands (the valence band 
of the electron-doped layer and the conduction band of 
the hole-doped layer) do not intersect the Fermi energy. 



Calculations described in this paper are based on a 
two-band model in which we retain only bands that have 
a Fermi surface; the role of remote bands is discussed 
for both topological insulator and graphene double-layer 
cases in Sec. |VI| In order to apply the two-band model 
consistently we impose an ultraviolet momentum cut-off 
at k = 2kp. 



We allow for the possibility that the bare Coulomb in- 
teraction strength is reduced by a finite fraction e due to 
extrinsic dielectric screening from the surrounding ma- 
terial. (A more elaborate model is required to account 
for the screening from the bulk of the TI as we discuss 
later.) Since both the band and interaction energies scale 
with distance L , this model has only three parameters 
when all energies are expressed in terms of the Fermi en- 
ergy ep : i) the Coulomb interaction strength parameter 
a = e 2 /etiVD ii) the Dirac band degeneracy N and hi) the 
dimensionless layer separation parameter kpd. Our main 
goal is to estimate how the spontaneous coherence gap 
at zero temperature depends on these three parameters. 



The starting point of any many-body theory is a 
derivation of the two-particle matrix elements of the 
electron-electron interaction between pairs of band eigen- 
states. We allow for the presence of a non-zero bare single 
particle tunneling amplitude, to, which will induce some 
degree of interlayer coherence even when the interactions 
are absent. Indeed, the interlayer tunneling transport 
anomalies on which BiSFET operation depends require 
at least a modest value for tJ3l. We will be interested in 
the many-body-enhanced tunneling amplitude, t, which 
includes a contribution from the electron self-energy that 
is off-diagonal in layer index. Spontaneous coherence oc- 
curs when t remains finite for to — > 0. Including to, the 
many-electron Hamiltonian is 
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where {pxkai^Xka} ar e the creation and annihila- tion operators for Dirac band states with definite two- 



dimensional momentum k and layer index X, while 
{ c k a > c ka} and {v ka ,Vk a } are creation and annihila- 
tion operators for the respective conduction and valence 
bands in the presence of the interlayer hybridization am- 
plitude t . In this equation E k = y/v 2 D {k — kp) 2 + £q is 
the bare band energy magnitude, 4>i is the orientation 
angle for momentum fcj, and kp is the common Fermi 
momentum of conduction and valence bands in the ab- 
sence of hybridization. The Greek indices on the creation 
and annihilation operators range over the number N of 
Dirac fermion flavors. When the Dirac model is viewed as 
a continuum model of graphene, we should choose A = 4 
to account for two valleys and two spin states, whereas 
for the surface states of topological insulators like Bi 2 Se 3 
we have N = M 

The linear dependence of band energy on momentum 
for to — in Hamiltonian (IT]) is a hallmark of the Dirac 
physics of graphene and topological insulators. The in- 
teraction term describes pairs of electrons which are scat- 
tered by their mutual interactions. In this term V x x' (?) 
is the two-dimensional Fourier transform of e 2 /r which 
equals Vs(q) = 27re 2 j q for electrons in the same layer and 
Vd(<?) = exp(— qd)Vs(q) for electrons in different layers. 



The factors of cos((0i— 4>a)/2) and cos(((j>2— 4>3}/2) which 
appear in this equation are well known in the graphene 
literature, and are due to the dependence of the phase dif- 
ference between sublattice components of the band states 
on the direction of momentum. 



III. BILAYER SCREENING 

In a bilayer, Coulomb interactions in one layer induce 
a charge response in the same layer and in the oppo- 
site layer. We use the random phase approximation for 
screening in which electrons respond like non-interacting 
particles to the sum of the external potentials and the 
mean-field Hartree potentials from the charge densities 
that they induce. The dynamic non-interacting response 
functions can be calculated by using first order time- 
dependent perturbation theory. When the responding 
electrons are coherent combinations of the two layers 
we find that the imaginary frequency same-layer and 
different-layer single-particle density response functions 



IL s (q,iu) = -N 
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v D {k - k F ), E k = ^Jv 2 D {k~k F ) 2 +t 2 , 



and t is understood as the full interaction enhanced hy- 
bridization amplitude. In solving the gap equations dis- 
cussed in the next section, we will find it convenient to 
take advantage of the smooth frequency dependence of 
the polarization functions along the imaginary axis. 

For an equal carrier density system, the screening of 
the sum V+ = Vs + Vd (even screening channel) and 
difference V- = Vs — Vd (odd screening channel) of the 
intralayer and interlayer interactions decouple. An el- 
ementary calculation yields simple geometric series ex- 
pressions for the even and odd screened interactions: 



V£ c = 



V± 



i-v ± {u s ±n D )' 



(4) 



Recalling that E 2 = ^ 2 + t 2 , it is easy to show that 
lis + IId vanishes for q — > 0. There is no even channel 
polarization and no even channel screening in this limit 
because conduction and valence quasiparticle bands at 
the same k are orthogonal. More fundamentally, the total 
density response of a gapped system to a uniform poten- 



tial shift is zero. For the q,ui — > limit, the odd channel 
perturbation on the other hand, corresponds simply to a 
shift in Fermi wavevector. It follows that the odd chan- 
nel polarization approaches the thermodynamic density 
of states in this limit, II5 — Ud — > —Nuq up to small 
corrections 0(t/eF) 2 , and it is not substantially altered 
by interlayer coherence (vq = kp /{2ttvd)). For q, ui —$■ 0, 
Vd —> Vs so that there is no odd channel potential, and 
screening is therefore absent as illustrated in Fig. 2. For 
larger q, Vs > Vd and the screened interlayer interaction 
can therefore sometimes be stronger than the bare inter- 
layer interaction. After a bit of simple algebra, we find 
the following general expression for the random phase 
approximation screened interlayer interaction: 
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The screened interlayer interaction is invariably larger 
than the value l/(2JWo) obtained when the polariza- 
tion functions are replaced by their zero frequency nor- 
mal state values. We therefore expect to find stronger 




in screening, than in the estimates of Ref. [TU] in which 
a static screening approximation for the polarizablizities 
in the normal phase was employed. 

IV. CONSTANT GAP APPROXIMATION 
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FIG. 2: Ratio of screened interlayer potential to bare inter- 
layer potential as a function frequency w along the imaginary 
axis. In these plots v is a non-linear dimensionless frequency 
defined by ui/tF = *V(1 — v) so that oj — €f <=> v = 1/2 and 
u = oo 4*- v — 1. Screening power increases with the prod- 
uct of the interaction strength a and the number of Dirac 
flavors N in each layer. The value Na = 8 corresponds to 
graphene surrounded by low e ~ 1 dielectric material. Screen- 
ing weakens at larger frequencies and larger gaps, vanishing 
in the high frequency limit as illustrated in panel (a) where 
the curves correspond to gaps A/ef = {0.1, 0.2, ..., 1.0} (The 
inset shows the zero frequency value of this ratio as function 
of the gap.). The flavor number N enhances screening as il- 
lustrated in panel (b). In isolation this consideration would 
favor enhanced coherence in the N — 1 topological insula- 
tor thin film case compared to the N — 4 graphene double- 
layer case. Screening vanishes for q,cu — > as illustrated 
in (c). The screened interaction can sometimes be stronger 
than the unscreened interaction as illustrated in panel (d), 
where the curves correspond to dimensionless layer separa- 
tions kpd — {0,0.4, 0.8, ..., 2.0}, because of the interplay of 
even channel and odd channel screening explained in the text. 



interlayer coherence in our calculations, which account 
for the retarded (frequency dependent) character of the 
Dirac gas polarization and include the effects of the gap 



In many-body theorjflSESl the influence of electron- 
electron interactions on electronic structure is captured 
by the self-energy. The self-energy can be thought of as 
an interaction correction to the quasiparticle Hamilto- 
nian, like the exchange and correlation corrections that 
appear in density-functional theory, except that it varies 
with the momentum and energy of the electron- wave un- 
der consideration. Our analysis of interlayer hybridiza- 
tion is based on the simplest sensible approximation for 
the self-energy in an interacting electron system. The 
approximation accounts for screening, which plays an es- 
sential role in any interacting electron problem in metals 
because it converts the long-range e 2 /r Coulomb inter- 
action into a shorter range interaction that has a finite 
value when integrated over space. It also accounts for 
the time-delay of screening responses, retardation, which 
makes screening less effective when inter-electron poten- 
tials fluctuate dynamically. We approximate the self- 
energy by the hrst order term in its perturbation theory 
expansion in terms of dynamically screened Coulomb in- 
teractions. This approximation is commonly employed in 
ab initio electronic structure theor y calc ulations where it 
is known as the GW approximatiorP^LLSl i n the graphene 
case this approximation has been employed to study a 
variety of physical properties both in the context of the 
Dirac equation electron gas niodePSl and in ab initio^* 
calculations and appears to be accurate where compar- 
isons with experiment are possible. One of the great 
successes of GW theory more broadly is its explanation 
for the property that band gaps in semiconductors are al- 
most always larger than predicted by density-functional- 
theory 1718 . This property is closely related to the en- 
hanced interlayer coherence on which we focus, which 
may be viewed as simply another example of this well 
known tendency. The interaction enhanced gaps depend 
essentially on the non-locality of screened exchange inter- 
actions which is neglected by the local density approxi- 
mation and in most exchange-correlation approximations 
that are commonly used in ab initio density-functional- 
theory calculations. 

The recipe for constructing the GW approximation for 
the self-energy from the Hamiltonian is well-established 
and leads in the case of a double-layer 2DEG Dirac model 
to a self-energy that is wavevector and frequency depen- 
dent and is a matrix in layer index: 



j x,x 



(«,»«) = 



d 2 k 1 



cos 



- <i>k+q 



du 

2^ 



G x ,x>(q + k, i(fi + w)) V£ x , (fc, iw), 



(6) 



where the Green's function is related to the self-energy 
through the Dyson equation: G — (Gq 1 — E) _1 . We have 
chosen to take advantage of the analytic properties of the 
dependence of screened interaction on complex frequency 
to evaluate this quantity along the imaginary axis, where 
its frequency-dependence is smoother than along the real 
one, and thus easier to treat numerically. 

The particle-hole and layer symmetry of our model 
leads to a self-energy matrix of the form 



Using the explicit form of the bare Green's function with 
interlayer tunneling amplitude to, 



G 1 (q,iuj) 



iu - £, q -t 
-t iuj + £,, 



(8) 



T,(q,iu>) 



_ (K{q,iuj) A(q,iuj) 
A(q,iuj) -K*(q,iu) 



we find that the off-diagonal component of the self-energy 
(7) satisfies the following equation: 



A(g,iu/) = 



d 2 k 1 + cos( 



{2-Kf 



duj to + A 

2^ (t + A) 2 + {u - ImKf + (£ fc + ReK)- 



V^{q-k,i{uji-uj)). (9) 



Throughout this study we neglect the interaction cor- 
rection to the diagonal self-energy K(q 7 iuj). Instead, we 
focus on the off-diagonal anomalous self-energy A(q,ioj), 
whose value determines the excitonic gap. For gaps 
much smaller than the Fermi energy the integral kernel 
in Eq. fcty is so highly peaked at the Fermi surface, that 
any trial gap function which only deviates from the ac- 
tual one at large momenta and frequencies will not af- 
fect substantially the self-consistent determination of the 
value of the gap at uj = and q = kp. In this section 
we discuss the simplest self-consistent procedure which 
exploits this property, in which A(q,iu>) is replaced by 
its value at uj — and q — kp, i.e. by its value at the 



center of the wavevector-energy region in which coher- 
ence is established. We will refer to this approximation 
as the constant gap approximation. 

The main merit of this approximation is that it allows a 
transparent discussion of the dependence of the coherence 
gap on the three dimensionless quantities, {a, kpd, TV}, 
that parametrize the problem. The constant gap approx- 
imation is relaxed in the next section, in which the co- 
herence gap is calculated more accurately by including 
its momentum dependence in the integral equation. 

In this approximation the interaction-enhanced inter- 
layer tunneling amplitude is t = to + A, where 



A = 



d 2 k 1 + cos(</> 9 — cf)k) 
(2tt) 2 2 



dw 



+t 2 +a 



r 



2 V^{\k-k F q\,iuj). 



(10) 



Here q is the direction of the wavevector on the Fermi 
surface for which we evaluate A; because of the isotropy 
of the problem A does not depend on q. 

Strictly speaking, bilayer exciton condensation is char- 
acterized by the spontaneous appearance of interlayer co- 
herence in the absence of any interlayer hybridization. 



Because single-particle processes that support interlayer 
hybridization are always present, strict exciton conden- 
sation is always a theoretical concept, but one which is 
still useful in considering real physical properties. 

The constant gap approximation is conveniently dis- 
cussed in terms of the quantity 



A f d 2 k l + cos(^-0 fc ) [°°dw 1 

t J {2ixy 2 J TT LO Z + t Z + ££ 



(11) 
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Coherence is possible in the absence of a single-particle contribution when /(t) = 1. This equation always has 
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FIG. 3: f(t) = A/t vs. t/eF on a logarithmic scale for N — 4 
and layer separation d — > at a series of values of a. The 
linear dependence of / on —Logio(t/ep) at very small gaps 
is a result of the nesting of electron and hole Fermi surfaces 
as explained in the text. The size of the spontaneous gap is 
determined by solving f(t) — 1. This figure illustrates why 
the spontaneous gap can be ~ ef at strong interactions and 
orders of magnitude smaller at weak interactions. 



a solution since f(t) — » for very large t, but diverges 
logarithmically for t — > 0. In Fig. [3] we plot f(t) vs t for 
N = 4 for a series of a values. 




Equation pi for the spontaneous gap is identical to f IG ' 4: a) Constant gap approximation spontaneous gap as a 



the equation for the superconducting gap in the BCS the- 
ory of superconductuctivity, except that in the Coulomb 
interaction case we must account for dynamic screening. 
The logarithmic divergence at small t can be traced to 
the nesting between conduction and valence band Fermi 
surfaces in bilayers with overall neutrality. Because the 
numerical coefficient in front of the logarithmic diver- 
gence is strongly reduced by screening^, /(£) reaches 1 
only for extremely small values of t when interactions are 
weak. 

The role of screening changes dramatically when inter- 
actions are strong. In Fig. [4] we have plotted the sponta- 
neous gap as a function of the dimensionless interaction 
parameter a = e 2 /evu in the limit of vanishing interlayer 
separation kpd — > 0, which is the most favorable limit for 
excitonic condensation. In graphene a steep increase in 
the gap occurs near a ~ 1.5, where it rises to become 
of the order of the Fermi energy. This result underlines 
the importance of the dielectric medium surrounding the 
double-layer having as low a dielectric constant, e, as 
possible. 

Self-consistent screening of Coulomb interaction in the 
presence of the excitonic condensate leads to the possi- 
bility of a first order quantum phase transition between 
states with very different values of the excitonic gap as 
the interlayer separation is changed. This originates from 
the tendency of the condensate to reduce screening, thus 
the screened interaction, V^ c , in Eq. (10) is an increas- 
ing function of A. At strong interactions this can lead to 



function of effective fine structure constant a — e /(svd), for 
(left to right) unscreened interactions and N = {1,2,3,4} 
degenerate Dirac cones (TV = 1 for Bi2Se3, TV = 4 for 
graphene). b) Comparison between several approaches to 
screening. When the screening is treated by using static 
screening an extremely small gap is predicted. If the reduc- 
tion in screening in the coherent state is neglected the sudden 
rise at a ~ 1.5 is absent (for the dashed line we fixed the gap 
to be 10~ 4 £f inside the polarization functions, thus virtually 
computing them in the gapless phase.). The full dynamical 
and gapped screening result approaches the result obtained 
with bare interactions at strong coupling but differs from it 
by several orders of magnitude at weak coupling (the crosses 
in the dynamic gapped screening case were obtained by ex- 
trapolating the / function to very small gaps see Fig. 131. 



non-monotonic behavior of the right hand side of Eq. ( 11 ) 
as function of A, as illsutrated in Fig. [5j 

When plotted as a function of interlayer distance for 
fixed a, the gap is a multivalued function for a window of 
distances. The branch on which it appears to increase as 
a function of distance at fixed densities and a is unphysi- 
cal, and the other two branches correspond to metastable 
states. Figure [6] illustrates this behavior for graphene 
with €f = 0.25eV and a series of a > 1.6. Thus at fixed 
a and large interlayer distances d 3> 1/fcjr the gap is sev- 
eral orders of magnitude smaller than ep, this tendency 
persists down to zero interlayer separation for small fine 
structure constants, but for a > 1.5, a first order phase 




FIG. 5: (color online) Evolution of f(t) with increasing in- 
terlayer distance. In these curves a = 1.8 and k F d = 
{0.04, 0.06, ...,0.12}. Metastability exists for k F d cr2 > k F d > 
k F d cr \, with k F d cr i ~ 0.06 and k F d cr 2 ~ 0.11. 



transition to a coherent state with a gap on the order 
of ep is encountered at a critical distance d cr < 1/kp. 
Metastability of the two states can persist down to a sec- 
ond critical value of the interlayer distance, after which 
only the state with strong interlayer coherence is a self- 
consistent solution of the gap equation. 

The same qualitative behavior is found for TIs as illus- 
trated in Fig. [7J although the first order phase transition 
becomes less pronounced at small N. Figure [7] also il- 
lustrates that the ratio A/ep is a decreasing function of 
kpd, this implies that for a given device, with a fixed 
interlayer distance and effective fine structure constant, 
the excitonic gap is not expected to be a monotonic func- 
tion of e^, but would increase at small ep, be peaked for 
a value ep < vp>/d, and decrease at larger values. 

The existence of a first order phase transition as a 
function of distance at zero temperature suggests the ex- 
istence of a first order transition as a function of tem- 
perature at fixed {a, kpd}, and thus a calculation where 
a continuous phase transition to the excitonic phase is 
assumed and the finite temperature linear gap equation 
is solved with screening computed in the normal phase 
would tend to underestimate the critical temperature for 
this transition. 

Spontaneous coherence in topological insulators is fa- 
vored by their reduced number of degenerate Dirac 
conea^Sl. On the other hand the 3D bulk of the mate- 
rial is expected to substantially screen the interactions 
between the surface Dirac electrons, strongly reducing 
their effective fine structure constant from its nominal 
value in vacuum ao = e 2 /vr>. In addition, the maximum 
possible value of the Fermi energy is the energetic sep- 
aration between the Dirac point and the closest band, 
which is ~ O.leV in current topological insulators. Even 
if large excitonic gaps were realized at zero temperature, 
the Kosterlitz-Thouless transition would limit the criti- 
cal temperature to ~ OTe^ ~ O.OleV^, this alone poses 
a major limitation to achieve ambient temperature inter- 
surface coherence in currently known topological insula- 
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FIG. 6: (color online) Spontaneous gap in graphene as a 
function of interlayer separation (with e F — 250 meV) for 
a = {1.6,1.8,2.0,2.2,2.4} (for graphene suspended in vac- 
uum a ~ 2.2, for an Si02 substrate a ~ 1.3). As discussed 
in the text there can be three self-consistent solutions at a 
given interlayer separation, the dashed line depicts the un- 
physical solution, the solid blue line depicts the physical so- 
lution where the gap reaches values on the order of the Fermi 
energy at small interlayer distances, and the red line schemat- 
ically depicts the solution where the gap is much smaller than 
the Fermi energy at large interlayer separations. The panels 
below present the schematic behavior of the free energy as a 
function of the gap, illustrating the first order quantum phase 
transition. 



tors. 

To assess the prospects of these materials more realisti- 
cally, we consider more closely one of the most promising 
candidates, namely Bi 2 Se^H It has a Fermi velocity of 
about 5 x 10 5 m/s which leads to an effective fine structure 
constant «o ~ 4.4 in the absence of dielectric screening. 
Unfortunately, Bi2Se3 has a large dielectric constant. To 
account for its influence we model the bulk of the topo- 
logical insulator as a uniform dielectric and assume that 
the thin film is surrounded by an extrinsic embedding di- 
electric. Then, the bare intra and interlayer interactions 
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FIG. 7: Spontaneous gap for a N = 1 topological insula- 
tor like Bi2Se3, as a function of the dimensionless interlayer 
separation k F d for a = {0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8}. In the 
absence of dielectric screening a ~ 4.4 for Bi2Se3 because its 
Dirac velocity is smaller than that of graphene. The interac- 
tion strength is reduced even in vacuum due to screening by 
the bulk of the topological insulator. The lower panels are 
for N = {2,3}. The first order phase transition is present at 
strong interactions but it is less pronounced as N decreases. 



the topological insulator (cti) to that of the surround- 
ing medium (eg). For very thin topological insulator 
films, the Dirac dispersion is modified by single-particle 
hybridization between top and bottom surfaces^. This 
limits the smallest thicknesses of thin films for which we 
can hope to observe interaction induced coherence effects 
to a few quintuple layers (Q£), d > SQL ~ 3nm. Con- 
stant gap approximation results obtained for the fixed 
ratio A/cp = 0.1 and d = 3nm are shown in Fig. [81 by 
fixing the ratio to be an order of magnitude smaller than 
the Fermi energy we have considered the limitation that 
the Kosterlitz-Thouless transition poses in order to pro- 
vide an estimate for the highest transition temperatures 
achievable in this material. 

From the results of this section we can say that the 
prospect for observing interlayer coherence in graphene 
and topological insulators is optimistic, provided that 
small interlayer separations can be achieved and extrin- 



sic substrate screening minimized. The prospects of 
approaching room temperature are more favorable for 
graphene than for TI thin films in spite of its larger num- 
ber of screening channels, mainly because of the possibil- 
ity of achieving larger Fermi energies but also because the 
TI surface state layers are separated by the TI bulk which 
screens interactions strongly. The energetically remote tt- 
bands will further favor i nterlayer c oherence in graphene 
as pointed out previously * 23 * 24 * 26 * 2 -^ and discussed further 
in Sec. IVTl 
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FIG. 8: Interaction induced gap in a thin film topological 
insulator, with parameters suitable for Bi2Se3. These results 
correspond to film thickness of do = 3nm ~ 3QL and the 
fixed ratio A/ef =0.1. For larger film thicknessess the self- 
consistent A would be modified by a factor of do/d. The 
dashed line indicates the room temperature energy scale. 



V. MOMENTUM DEPENDENT GAP 

A simple way to include the momentum dependence 
of the gap is to cast Eq. ([9| as an eigenvalue problem. 
To do so, we replace the gap by its value at fl = and 
q = kp in the energy denominator and in the screening 
function while retaining full momentum dependence in 
the numerator of the anomalous Green's function. This 
procedure is justified if we assume a smooth wavevector 
dependence of the gap near the Fermi surface over a scale 
A/vp, and a weak frequency dependence over a scale ep, 
below we illustrate that these two conditions are met. 
With these simplifications the gap equation becomes an 
eigenvalue problem 



A(q) = 



d 2 k 1 + cos( 



(27T) 



did 



A(fc) 



7T w 2 + a 2 + e 



V^{\k-k F qllUj). 



(14) 



This equation is solved when an eigenvalue of the integral 
operator is equal to 1. By demanding that the solution 



is a fixed point of the integral equation, we conclude that 
the physically relevant solution corresponds to the largest 



9 




FIG. 9: a) Value of the momentum dependent excitonic gap 
at q — fe_F as a function of the effective fine structure constant, 
determined from the eigenvalue problem described in the text. 
This figure corresponds to the limit d — >• and the curves are 
for different Dirac cone degeneracies N = {1,2,3,4}. Notice 
that in graphene (N — 4), the gap is not single valued as a 
function of a indicating that a discontinuous phase transition 
can occur even for fc^d = when the momentum dependence 
of the gap is taken into account (compare with the results 
within the constant gap approximation of Fig. [4n. The inset 
shows the largest eigenvalue of the integral operator for A = 
O.IeF as a function of a for N = 4, self-consistency is achieved 
at an a for which A m — 1. b) Momentum dependence of the 
gap corresponding to the same parameters used in the inset 
of a), c) Frequency dependence of the gap as a function of 
the non-linear frequency scale, cj/ef = vj{\ — u), estimated 
by replacing the momentum dependent of gap of panel (b) 
in the L.H.S. of Eq. [9] The smooth frequency-momentum 
variation across the fermi surface over a scale A/vd validates 
the consistency of the approach. 



eigenvalue Am which must equal 1. By computing this 
eigenvalue as a function of the value of the gap at the 



Fermi surface, which enters the integral kernel through 
the energy denominator and the screened potential, we 
can solve for Am (A) — 1 to find the self-consistent value 
of the gap at the Fermi surface. In this approach Am 
therefore plays a role similar to that played by f(t) in 
the constant gap approximation. The eigenvector corre- 
sponding to the largest eigenvalue, normalized to have a 
value A at q = kp, is the momentum dependent gap func- 
tion. Compared with the constant gap approximation we 
find that accounting for the momentum dependence of 
the gap reduces its value at the Fermi surface, although 
the qualitative trends towards coherence obtained within 
the constant gap approximation are the same, so that the 
conclusions of the previous section largely stand. 

Note that because we include A in the energy denom- 
inator and in screening, we are not solving a linearized 
gap equation like the one used in BCS theory to estimate 
the critical temperature. Non-linearity is encoded in the 
self consistency of the largest eigenvalue. This approach 
is justified when the imaginary frequency dependence of 
the gap is weak and gaps are smoothly varying over a 
width 5q ~ A/vd around the Fermi surface, since this 
is the momentum region that contributes to the integral 
kernel. Both assumptions are justified as illustrated in 
Figs. |9}>) and [§;). 



VI. CONTRIBUTIONS FROM REMOTE BANDS 

So far we have been dealing with a simplified two-band 
model to show that a self-consistent inclusion of the ex- 
citonic effects into the screening properties of a bilayer 
system considerably enhances its tendency towards in- 
terlayer coherence. Further, we used a hard momentum 
cut-off at 2kp in the calculations. In this section we dis- 
cuss what would be the effect of including the momenta 
beyond our Ikp cutoff, and the remote Dirac bands which 
are not crossed by the chemical potential. Generally two 
competing factors appear, the enhancement of screening 
from the addition of available states for single particle 
excitations, and the increment of available phase space 
for coherence with a tendency to raise the gap as we will 
describe in this section. 

Consider the full four band model equations with band 
diagonal pairing (see e.g. Eq. (3.1) of refP^j): 



A s (q,u) = 



Z^ 



d 2 q'du' 
(2tt) 3 



Vtftq-q'^-Lo')) 



ss cos 



uj ■ 



A 



A, 



■a 



(15) 



r 



where s = {+,—}, £ s = vp>(q + sfcp) and the ' accent 
on A indicates that its arguments are u)',q'. The point 
we wish to make can be recognized by considering the 
limit kp — > of this equation for which A + = A_ = A. 



The gap equation for spontaneous coherence in this limit 
is identical to the spin-density wave gap equatiorPSl of a 
single neutral graphene layer, except that in the present 
case the interaction is an interlayer rather than an in- 



tralayer interaction: 
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A(g,w) = 



dPq'duj' , , 

, 9 ^ Vg c (g - q', iw - «')) — 

(27TJ J w 



A 



2 '9 ' 



(16) 
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In the single-layer case it is known that beyond a crit- 
ical interaction strength a gap will open due to coher- 
ence between conduction and valence bands that marks 
the onset ol a transition to a gapped state with sub- 
lattice spin-density wave ordering^. (This ordered state 
is sometimes referred to as an excitonic state, although 
this terminology is not appropriate since the number of 
states in the valence and conduction bands of a single 
layer are not separately conserved, even approximately.) 
Since the neutral system lacks an energy scale, the size of 
such gap is determined by ultraviolet scales beyond the 
Dirac model. Just as a spin-density wave state is a real- 
istic possibility in neutral single-layer graphene, an inter- 
layer coherent state is a realistic possibility in graphene 
double-layers even in the absence of carriers. States that 
lie above the 2kp ultraviolet cut-off we have employed so 
far can induce coherence on their own if interlayer inter- 
actions are strong enough, so they must have an influence 
even at more moderate interaction strengths. 

To illustrate their influence, we separate the integral 
in Eq. ( 15 1 into two contributions, one for momenta run- 



ning up to 2kp and one for momenta running from 2kp 
up to a new UV momentum p w on the Brillouin-zone 
scale. For simplicity we assume that A + = A_ = A; 



it is easy to show that this becomes exact for the case 
of a pseudo-potential Vp c (q — q' ,i(u> — u/)) which is in- 
dependent of the angle between q and q' . (Note that 
A+(<7 = 0,iui) = A_(g = 0,iu>) even without this as- 
sumption, and that the coherence of remote bands con- 
sequently does not necessarily open up a gap at the Dirac 
points of each layer.) 

We approximate the large momentum contribution to 
the right-hand-side of the gap equation by assuming that 
the gap varies as 1/q at large q because of the large q 
behavior of the screened Coulomb interaction, writing 
A(q,itu) ~ Ap(kp/q), where Ap is the value of the gap 
at the Fermi surface. Calling SAp the contribution to 
the gap at the Fermi surface coming from UV momenta, 
we can write for p w ^> kp, 



SA f 



2k, 



(27T) 
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The frequency integral can be approximated by view- 
ing it as averaging the screened potential with a non- 
normalized Cauchy-Lorentz distribution of width ~ Vp>q'- 



duj V£ c (q,iu) 1 

n u 2 + (v D q) 2 2voq 



(Vl, c (q,ioj))« 



Vi, c (q,Q(vDQ)) 

2v D q 
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In the limit d — > 0, it follows from the behavior 
of the polarization function in neutral graphene^ that 
Vj} C (q,0(vDq)) — > (27re 2 )/g, where e 2 is a function of 
a, N which measures the effective strength of the inter- 
action at high energies. In this way we obtain, 



SAp 



a* 



A, 



a* 
a c 



(19) 



where a* = e 2 /vr> and a c is the value of a* for which 
coherence would occur without carriers. Although this 
line of argumentation is qualitative, two points seem un- 
avoidable: a) that the contribution to the right-hand 
side of the gap equation from UV momenta is propor- 



tional to the value of the gap near the Fermi surface, 
and b) that the factor of proportionality is smaller than 
one unless interlayer coherence would occur even without 
carriers. The absence of a spontaneous gap in neutral 
single layer graphene^ suggests that spontaneous coher- 
ence gaps should also be absent in neutral double-layers 
and therefore that the factor of proportionality should be 
smaller than one. The UV remote band contribution to 
spontaneous coherence at finite carrier density can nev- 
ertheless be substantial. By adding the large momentum 
contribution to the one coming from momenta p < 2kp 
we can rewrite the gap equation in the form 
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The UV contribution can be captured by multiplying 
the screened interaction in the gap equations we have 
used in previous sections by an enhancement factor S = 
(1 — a±/a c )~ l > 1. This analysis could be thought of as 
a renormalization transformation that maps the original 
model defined up to momentum scales of p w into a model 
of momenta up to 2kp while keeping the Green's func- 
tions fixed at the expense of rescaling the interactions 31 . 
Of course the exact transformation would also alter the 
frequency/momentum dependence of the interactions. 

It follows from Eq. (t) that SV£, c (q,iui;a, N,d) = 
Vf) C (q,iuj; Sa,N/S,d). The effect of the interaction en- 
hancement is therefore two-fold, to increase the effective 
fine structure constant a — > Sa and to reduce the effec- 
tive degeneracy of the Dirac cones N — > N/S. Both ef- 
fects support coherence and their influence could be quite 
substantial. These considerations motivate the two-band 
calculations in Fig. [4] for values of a and N that do not 
correspond directly to either TI's or graphene. 

In the right side of Eq. ( 20 1 the bands that cross the 



Fermi surface have a singular contribution in the limit 
Ap — >• which guarantees a solution. The remote band 
contribution does not on its own guarantee a solution 
but, when included, can only increase the estimate of 
the gap. This prope rty has already been emphasized in 
previous studiea^EH. 

We do note again that the polarization functions that 
enter Eq. ( 14 ) are those from the four-band model, 



whereas we have employed two-band polarization func- 
tions in our two-band model calculations. The polar- 
ization functions of the two models coincide for mo- 
mentum transfers q < kp, but at higher momenta the 
four-band model polarizations produce stronger screen- 
ing. Nevertheless, the contribution of such momenta to 
the two-band gap equation is reduced by the backscatter- 
ing suppression factor of Dirac models and by the decay 
of Coulomb interactions at larger momentum transfers. 

We believe our approach captures the qualitative 
trends in the cxcitonic condensation in double-layer sys- 
tems, although precise quantitative predictions would re- 
quire a more elaborate numerical approach which incor- 
porates the competition of the effects described in this 
section. 



VII. SUMMARY AND DISCUSSION 

This article addresses the possibility of achieving spon- 
taneous coherence between the layers of a graphene 
double-layer system, or between the top and bottom sur- 
faces of a topological insulator thin film at temperatures 



needed for its robust experimental observation, or even 
for practical applications. 

The task is accomplished by solving for the T = 
energy gap A using a GW dynamic screening approxi- 
mation for interlayer interactions. The accuracy of this 
approximation in single-layer graphene systems is sup- 
ported by detailed comparisons^ with ARPES spectra. 
Although we believe that the approximation we employ 
is valuable for the purpose of identifying the most favor- 
able circumstances for interlayer coherence, it is not fully 
rigorous. For example, it neglects the rearrangement of 
spectral weight due to quasiparticle renormalization fac- 
tors that is commonly included in strong-coupling theo- 
ries of superconductivity. More seriously, in our opinion, 
it does not adequately capture the competition between 
coherent states, which lower interlayer interaction ener- 
gies, and other strongly correlated states which minimize 
intralayer interaction energies. 

The two double-layer systems that we focused on are 
distinguished by the number N of Dirac fermion flavors, 
by their Dirac velocities vd , and by the dielectric screen- 
ing environments that surround the two-dimensional itin- 
erant carriers. Spontaneous coherence in these sys- 
tems can be promoted by using gates to induce two- 
dimensional electron and hole gases with carrier densi- 
ties that are equal in magnitude, but opposite in sign in 
the two layers. Because of the particle-hole symmetry 
of two-dimensional Dirac systems, this circumstance will 
guarantee accurate nesting between the electron-like and 
hole-like Fermi surfaces. Under these circumstances, it 
is generally thought that repulsive Coulomb interactions 
between electrons in different layers are guaranteed to 
induce spontaneous coherence at sufficiently low temper- 
atures. For small ratios of the excitonic gap, A, to the 
carrier Fermi energy ejr, mean-field theory is expected to 
be accurate at finite-temperatures. Recalling the famil- 
iar result from the BCS theory of superconductivity^, 
this implies a critical temperature ksT c ~ 0.5A. At 
larger values of A the phase transition, which is of the 
Kosterlitz-Thouless type, will have a critical temperature 
that saturates. Because the Kosterlitz-Thouless temper- 
ature cannot be larger tharP ~ O.lep, we are interested 
in identifying circumstances under which A ~ O.Ief or 
larger. Under these circumstances, achieving high tem- 
perature coherence is simply a matter of increasing the 
carrier Fermi energy. 

When the retarded character of the screened inter- 
layer interaction is ignored by assuming the statically 
screened form of the Coulomb interaction to be valid 
up to frequencies of the order of the Fermi energy, it 
can be shown that the gap can, at best, reach scales of 
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A < ep exp(— 4iV)i2], where the flavor number N = 4 
for graphene and N = 2 for topologic al in sulator sur- 
face states. Since the largest achievable^ Fermi ener- 
gies are smaller than ~ leV in current graphene samples 
and smaller than ~ 0.1c V in current topological insu- 
lators, the largest achievable gaps would be ~ 10 -7 eV 
in graphene double-layers and ~ 10 _3 el^ in topologi- 
cal insulators if this approximation was accurate. Our 
main finding, which is summarized in Fig. HI is that a 
self-consistent description of dynamical screening in the 
gapped phase leads to zero temperature gaps which be- 
come of the order of the Fermi energy when the layers are 
close together and the dimensionless coupling constant 
which characterizes interactions strengths a = e 2 /eHvp, 
is large. 

The double-layer electron-hole pairing problem shares 
some features with the BCS model of electron-electron 
pairing in conventional superconductors. There are nev- 
ertheless important differences, one of the most promi- 
nent ones being the lack of an intermediate energy scale 
between the gap and the Fermi energy which in con- 
ventional superconductivity is provided by the Debye 
frequency. In the present case the screened interac- 
tion crosses over as a function of frequency from the 
statically screened values at zero frequency to the un- 
screened Coulomb potential at frequencies on the order of 
the momentum-transfer dependent plasmon energies^, 
which vanish at large wavelengths in two-dimensional 
electron systems. Since the momentum transfer for 
scattering at the Fermi surface can take values q e 
[0,2kp] 7 the small momentum transfer region is partic- 
ularly poorly described by a statically screened interac- 
tion. This effect is responsible for increasing the gap by 
orders of magnitude 11 compared to the non-zero but un- 
observably small values, obtained with pure static screen- 
ing!^. An additional enhancement is brought about by 
self-consistently accounting for the effect of the gap on 
screening. The gap suppresses screening only at small 
momenta q < /S./vp>, but this is where the Coulomb in- 
teraction is strongest. Therefore accounting for the de- 
crease of screening when the gap opens makes a sizable 
contribution to the right hand side of the gap equation. 
These different approximations are compared directly in 
Fig. [4] The final result is that the gap changes relatively 
suddenly, sometimes discontinuously, from a value that 
is orders of magnitude smaller than ep to a value larger 
than O-le^ at a critical interaction strength which is close 
to the values that can be achieved by minimizing dielec- 
tric screening. 

Self-consistent screening of Coulomb interaction in the 
presence of an excitonic gap has important consequences 
for the fate of the condensate as the interlayer distance is 
increased. In particular, we found that a first order phase 
transition from a large gap state to a small gap state oc- 
curs as interlayer separation increases. In the presence 
of disorder the smaller-gap phase would be hard to dis- 
tinguish from a normal phase, thus in experiments this 
transition should effectively look like a first order transi- 



tion between coherent (large-gap) and normal phases. 

The results in this work were obtained by consider- 
ing only those two-dimensional bands which have Fermi 
surfaces. In the case of topological insulator thin films 
there are many bulk-like bands at higher energies. It 
seems likely that the only important role of these bands 
is in providing additional screening which can be approx- 
imately captured by modeling the TI bulk as a uniform 
dielectric. 

The remote Dirac bands, that is the valence band of the 
n-type layer and the conduction band of the p-type layer, 
play a special role in establishing the interlayer coherence 
in graphene. Their supportive role, described in more de- 
tail in Sec. |VI[ can be understood from the following qual- 
itative argument. In the two-band model the self-energy 
responsible for the order of the condensed state induces 
coherence between the conduction band of the n-type 
layer and the valence band of the p-type layer. In the lan- 
guage of graphene's sublattice pseudospin (or real spin in 
the case of TIs) such a self-energy represents an interlayer 
coupling between antiparallel pscudospins. Because the 
remote band states in each layer differ from those of the 
bands with Fermi surfaces only by pseudospin reversal, a 
self-energy that couples one set of bands will also couple 
the other set. The two contributions to interlayer coher- 
ence are in phase at all wavevectors. The remote bands 
can have a large impact, even though they are removed 
energetically, because they are present over a very wide 
energy range. The pseudospin reversal relationship be- 
tween conduction and valence bands in graphene which 
leads to this property is paralleled by a real spin reversal 
relationship in topological insulators. The remote band 
effect therefore applies equally well to that case, but over 
a much more limited energy range. The quantitative im- 
pact of remote bands is difficult to estimate, however, 
as discussed in Sec. IVII If coherence does occur in the 
range of interaction strengths and layer separations that 
is achievable experimentally, the remote band effect could 
be essential. 

Our results indicate that reducing extrinsic screen- 
ing is crucial for the observation of interlayer coher- 
ence. Consider the case of hexagonal boron nitride (h- 
BN), an attractive material for building high quality 
graphene hetero-structures^l. h-BN has a dielectric con- 
stant e ~ 5 33 , thus, double-layer graphene fully embed- 
ded in h-BN would have an effective fine structure con- 
stant a ~ 0.4, and our estimates suggest that the gap 
will only rise to values A < 10~ 5 ep ~ 0.1K at ep ~ leV. 
On the other hand, for a heterostructure where h-BN 
is only used as a substrate and as a thin barrier while 
the top layer is surrounded by vacuum we have an ef- 
fective dielectric constant e ~ 3, for which a ~ 0.7 and 
A < 10~ 4 e^ ~ IK. Finally, a heterostructure where h- 
BN is used only as a thin barrier between the layers while 
the upper and lower halves are surrounded by vacuum 
should have an effective fine structure constant near that 
of graphene in vacuum a ~ 2.2, and our results indicate 
that in such case the zero temperature gap will become 
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of the order of the Fermi energy and interlayer coherence 
could survive to ambient temperatures. 

After completing this work we became aware of a re- 
cent study which reaches similar conclusions to the ones 
we have presented^. 
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